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1 Introduction 



The 1/iV, or topological, expansion in U(N) invariant matrix models is being considered with re- 
newed interest in the last years because of its various applications in topological string theories on 
special classes of Calabi-Yau geometries or J\f = 2 superconformal gauge theories [1, 2, 3]. The 
solution in the limit N — > oo is generically described by a complex curve. This is so also in the 
case when the symmetry group is a direct product of several unitary groups, as the ADE matrix chains 
[4, 5]. Since the complex curve determines the classical spectral density of the matrix variables, it 
is also called spectral curve. In all solved examples, the l/N expansion can be formulated entirely 
in terms of the complex curve. The fact that the spectral curve determines not only the classical 
limit N — > oo, but also the complete l/N expansion, is highly non-trivial and deserves to be better 
understood. 

Two approaches to the l/N expansion formulated in terms of the spectral curve were developed 
in the past ten years. The first approach, which will be refered to as CFT approach, is based on the 
conformal invariance, which is believed to be present in all U(N) invariant matrix systems. The 
conformal symmetry is obvious for the class of matrix models which, after diagonalization, reduce 
to Coulomb gases [6, 7, 8]. A prescription to evaluate the quasiclassical expansion in such matrix 
models using the toolbox of CFT was outlined by one of the authors in [8, 9]. The basic idea is that 
the conformal invariance is sufficient to construct the collective field to all orders in the quasiclassical 
expansion. The spectral curve appears as the classical value of a bosonic field, which lives on a Rie- 
mann surface associated with this spectral curve. The gaussian approximation gives the two leading 
terms of the l/N expansion, while the higher terms are obtained by inserting special local operators at 
the branch points of the Riemann surface. 2 The proposal of [8, 9] is a field-theoretical formulation of 
the 'method of moments' for performing the higher genus calculations, developped in [11], and gen- 
eralizes the CFT description of the non-critical string theories [12, 13, 14]. The case of a hyperelliptic 
Riemann surfaces was considered in details in [15]. 

Another approach to the l/N expansion, known as topological recursion, was developed more 
recently by B. Eynard and collaborators [16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26]. Whilst the CFT 
approach is related to field theory, the topological recursion, which we denote shortly by TR, is related 
to algebraic geometry. The topological recursion gives a very efficient algorithm for calculating the 
l/N expansion of the observables representing symplectic invariants on the complex curve. The 
recursion procedure overcomes some of the technical difficulties of the method of moments [11]. The 
basic observation here is that recurrence equation is most simply formulated in terms of residues at 
the branch points of the Riemann surface. The outcome of the topological recursion is an elegant 
graphical scheme, containing only trivalent vertices. 

Let us emphasize that the CFT and TR approaches can be formulated only in terms of the spectral 
curve, without any reference to the underlying matrix model. Therefore, instead of speaking of the 
l/N expansion, we will speak of the quasiclassical expansion (in h = l/N) associated with the 
spectral curve. Since any spectral curve can be considered as a solution of the universal Whitham 
hierarchy [27], CFT and TR give two methods to evaluate the quasiclassical expansion of an integrable 
hierarchy, given its dispersionless limit. 3 

2 The field-theoretical description of the branch points was pioneered by Al. Zamolodchikov [10]. 

3 Reconstruction of the quasiclassical expansion for the KP hierarchy, named ft-dependent KP hierarchy, was considered 
in [28, 29]. A serious problem in these works, which are based on the Lax formalism, is that the recursion relations are 
extremely complicated. The advantage of the CFT/TR method, based on the Virasoro symmetry, is that the solution is 
obtained in a simple and easy to manipulate form. We are not going to discuss here the uniqueness of the quasiclassical 
expansion. CFT/TR is believed to give the same quasiclassical expansion as the Lax formalism, but as far as we know there 
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The aim of this paper is to compare the two methods and the diagram techniques they lead to. 
Since the CFT method and the topological recursion represent two different techniques to resolve 
the Virasoro constraints, they are expected to lead to the same solution, although this is far from 
obvious for general spectral curve. On the other hand, the two diagram techniques are seemingly 
very different and have different geometrical interpretation in terms of world sheets. Therefore we 
think that it is important and instructive to find out the exact correspondence between the two forms 
of the quasiclassical expansion associated with a spectral curve. In this paper we show that the two 
diagram techniques lead to the same result and that moreover one of them is obtained by a partial 
resummation of the other. We conclude that the CFT method and the topological recursion are two 
different realizations of the same procedure, which we will refer to as CFT/TR. We underline that we 
are considering the case of the most general, not necessarily hyperelliptic, compact complex curve 4 . 

The paper is organized as follows. In Sect. 2 we collect some definitions and notations related to 
the spectral curve. In Sect. 3 we describe the CFT approach. We first remind the operator formalism 
for a gaussian field on a Riemann surface. Then we express the Virasoro constraints the operator 
formalism for a gaussian field on the spectral curve. We write the Virasoro constraints in terms of the 
mode expansion of the gaussian field near the branch points and construct the dressing operators to 
be inserted at the branch points in terms of these modes. Finally we derive the diagram technique for 
computing the higher genus free energies and correlation functions. In Sect. 4 we first give a brief 
summary of the topological recursion. Then we rewrite the equations of TR in terms of the series 
expansions of the symplectic invariants near the branch points. We start with the simplest case of a 
spectral curve with only one, simple, branch point. In this case the symplectic invariants are given by 
the observables of the Kontsevich model. We formulate the equations of TR in terms of a trivalent 
graph expansion. After that we consider the case of an arbitrary spectral curve. Sect. 5 we give the 
exact correspondence between the diagram techniques in the CFT and the TR approaches. 



2 Spectral curve and classical free energy 

By spectral curve we understand the triple 

£ = (X,x,y), (2.1) 

where S is a Riemann surface of genus G and the functions x and y are analytic in some open domain 
of E. The Riemann surface is assumed to have / punctures (/ > 1), where the differential ydx becomes 
singular. For later convenience we introduce a global parametrization z of the Riemann surface S. 
The complex variable z belongs either to the Riemann sphere (for genus 0) or to a quotient of the unit 
disk (for genus 1 or larger). In terms of the global parameter z we have 

x = x(z), y = y(z) , (2.2) 

where y and x are analytic functions of z. In the neighborhood of the punctures z = a a we chose 
local coordinates C(z,a) sucn that 

C ( aU=° (a = l,...,i). (2.3) 

is no general proof of that. 

4 From this perspective, the result of this paper goes beyond the unitary one matrix model. It applies for example ti the 
unitary two matrix models where the spectral curve is not hyperelliptic in general. 
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We can turn the spectral curve into a set of algebro-geometric data for Krichever's universal Whitham 
hierarchy [27]. We introduce the classical collective field as the one-form 



dS{z) = y(z)dx{z), (2.4) 

which is by definition holomorphic on the Riemann surface except at the punctures z = a a . The 
function S(x) is nothing but the potential for the universal Whitham hierarchy. The classical free 
energy associated with the spectral curve is defined as [27] 

^[S] = [ dSAdS. (2.5) 



The deformations of the spectral curve can be represented as commuting flows in the phase space 
corresponding to the "times" associated with the punctures. The Poisson brackets defined in the space 
of functions of the two variables x and y are defined as 

{/, g} = d x f d y g - d x g d y f (2.6) 

and correspond to the symplectic form uj = dx A dy. The classical field S can be thought of as the 
action of a Hamiltonian system depending on the coordinate x and the "times" associated with the 
expansion at the punctures [27]. 

The classical problem in the universal Whitham hierarchy is to determine the one-form (2.4) by its 
asymptotics at the punctures and the moduli associated with the non-contractible cycles. The moduli 
of the spectral curve at the punctures are given by the principal parts 

t a = Res y{z)dx{z) (2.7) 

z— >a a 

and the potentials 

V a (z) = Res y(z')dx(z')\n (l - . (2.8) 

z'^a a \ C(z',a) J 

Let {Ai, Bi}f =l be a canonical basis of cycles on E. The classical solution is completely determined 
by its asymptotics at the punctures and the 'filling fractions' 

v k = ^- j y{z)dx{z) (k = l,...,G) (2.9) 

A- 

associated with the A-cycles. A basis of holomorphic one-forms on the Riemann surface associated 
with the cycles Ak is given by the derivatives 

d 

iv k = -K—ydx. (2.10) 
dv k 

Their integrals along the A- and B -cycles are given by 

— <(> LOl=5 k l, — <L LOl=T kl , (2.11) 

2™ J Ak 2vn J Bk 

where r = {r^} is the period matrix of the Riemann surface. We will consider the filling fractions as 
a set of external parameters, which are part of the moduli of the spectral curve. 
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An important role in the construction of the quasiclassical expansion is played by the Bergman 
kernel for the spectral curve, B(z, z'). The Bergman kernel is denned uniquely by the following three 
conditions: 

- it is globally defined on the Riemann surface, including the punctures; 

- it behaves at small distances as 



dz dz* 

B(z,z') = —7T + regular function; 



(z - z<) 

- it has has vanishing integrals around the „4-cycles 



<p B(z,z') = 0, 



u 



1. 



,G). 



(2.12) 



(2.13) 



The Bergman kernel (2.12) depends on the spectral curve only through the 3G — 3 + 1 (for G > 2) 
complex moduli of the punctured Riemann surface S. It is important noticing that it does not depend 
on the moduli associated with the singular behavior of ydx near the punctures. 

Below we will assume that the classical problem is solved and the one-form (2.4) is already known. 
The free energy (2.5), or the classical action for the collective field, can be written as a sum of residues 
associated with the punctures and contour integrals of the one-form ydx [27, 30, 19]. For the sake of 
completeness we remind the general expression for J^ ) [S] , although we will not use it: 



^[S] = - 1 - 



G 

Res V a (z) y{z)dx(z) + V" t a fi a + v% <p y(z)dx(z 

^ z->a Q ' ^— ' J 

a a i=l „ 



(2.14) 



where 



d((z,a) 



Ha= / y{x)dx{x) - dV a (z) + t u 

'(z=a a ) V Q{z,a) 



o,aj 1 



(2.15) 



for an arbitrary base point o. 



3 The CFT approach 

The CFT approach is based on the assumption that the quantum collective field, which we denote by 
<1>, has the same analytical properties as the classical solution S(x). In other words, the field is 
invariant with respect to the conformal transformations of the spectral parameter x that preserve the 
I punctures x a = x(a a ) of the Riemann surface. 5 The outcome of the CFT approach is an universal 
formula for the quasiclassical expansion of the free energy 

JF[S, h]=Y^ tf 9 ~ 2 F {9) [S] + non-perturbative terms. (3.1) 

in terms of the classical solution S(x). 

5 The invariance with respect to conformal transformations which do not preserve the punctures determines the reaction 
of the system to a change of the external parameters. This symmetry is not relevant for the quasiclassical expansion. 
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3.1 Gaussian field on the spectral curve 

In the gaussian approximation, the classical field Q c \ is the expectation value of a gaussian field <£. 
We thus split the collective field into classical and quantum parts: 



The subleading term in the quasiclassical expansion for the free energy, -Fgauss = comes from the 
gaussian fluctuations around the classical solution. The corresponding factor in the partition function 
is given by the inverse power of the holomorphic piece of the determinant of the Laplace operator on 
the Riemann surface: 

-2gauss = —s ■ (3-3) 

detao 

The expression of the chiral determinant for a hyperelliptic curve has been given in [10] and for a 
general surface in [31, 32]. 

All correlation functions in the gaussian approximation can be expressed through the two-point 
function of a gaussian field on the Riemann surface, which is given by the Bergmann kernel, 

B(z,z') = <^ qu (z)^ qu (z')) s - (3-4) 

3.2 Operator formalism for the gaussian field on a Riemann surface 

Let zo be a base point on the Riemann surface £ and be a local coordinate in the neighborhood of 
zq such that £( zo ) = 0. We can associate with this variable a Hilbert space by expanding the gaussian 
field in Laurent series. We split the quantum field into a singular and a regular parts 

d$ qu (z) = d$ qn (zf s + d<S> qu (zf\ (3.5) 

having mode expansions 

n>0 n>l 

The amplitudes J n are assumed to satisfy canonical commutation relations 

[J m , Jn] = m5 m+nfi , m,n£Z. (3.7) 
This mode expansion defines a Hilbert space with left and right Fock vacua satisfying 6 

(0|J_ n = 0, J„,|0) = (n>l). (3.8) 
From the commutation relation 

[d*v{*r,d*zw]= f iz)d f z '\ 2 , (3.9) 



6 The amplitude Jo commutes with the rest and its action on the vacuum states can be specified separately. We can define 
a vacuum state with charge N by requiring that Jo |0) = AT|0), Alternatively one can introduce the zero mode <l>o. which is 
canonically conjugated to Jo, so that [Jo, $o] = 1- Then the vacuum state with charge iV is e^* |0). 
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one finds for the correlation function 

(0\d% u (z)d$ qa (z')\0) = r< 5 (2)d f (2,) 12 - (3.10) 

The left and right Fock vacua defined above are such that the gaussian field <I> qu has vanishing 
expectation value and its two-point function (3.10) gives the singular part of the Bergman kernel. We 
would like to deform the left and right vacuum states in such a way that the two-point function of <£ qu 
reproduces also the regular part of the Bergman kernel (2.12). In order to obtain the operator repre- 
sentation of the collective field <I>, eqn. (3.2), we should also introduce the expectation expectation 
value <fr c i. For that we first expand the classical field and the Bergman kernel in the local variable & z \, 

assuming that <E>™ 8 (zq) = 0, i.e. zq is a regular point on the Riemann surface: 



d$d(*) = lE^^y 1 ^)' (3 - n) 

n>l 

B M = l^-tl^ D B^nC^t^d^d^y (3.12) 
L? M ? ( 2 ) J n,m>l 

Then, as suggested in [33, 34, 35], we Bogolyubov transform the left Fock vacuum, 

<0| -+ <E| = (0| exp I I V ^JmJn + ^ E - J « I • ^ 

\ m,n>l n>l I 



In this way the Riemann surface E is represented by a state (E| belonging to the Fock space associated 
with the base point z$. For any correlation function of d$ qu (zj), with z\ belonging to a coordinate 
patch which contains the point zq, the unnormalized expectation value on the surface E is given by 

(d$qu(zi) ■ ■ ■ d% u (z n )) s = md<f> qu ( Zl )...d% u {z n )\0). (3.14) 

If we need to evaluate the correlation functions of operators that belong to different coordinate 
patches of E, we should generalize the operator representation for a Riemann surface with several 
extra punctures z\, . . . , z r . To each puncture we associate a pair of vacuum states as above and 
define the right and the left Fock vacua as the direct products of the Fock vacua associated with each 
puncture. Let C(z,s) t> e the local coordinate in the s-th coordinate patch. Then (3.13) generalizes to 



<0|=<g><0M| (S| = (0|exp[l £ ^fl^)jM + |^I^)j<.A 

s=l \ m,n>l n>l / 



(3.15) 



where 



m£ s) = Res y{z)dx{z) £ " 

Z — >Z a 

= Res Res B(z,z') £7^7?.. 
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3.3 Conformal invariance at the branch points 



The stress-energy tensor for the conformal transformations of the global parameter z is 

T(z)dz 2 = - lim \d<f>(z)d<f>(z') - B(z, z')] . (3.16) 

2 z'— YZ 

The gaussian field is by definition invariant under such transformations in the sense that the expecta- 
tion value ( T(z) ) is analytic everywhere except at the punctures z = a a where the classical solution 
has poles. 

The basic assumption of the CFT approach is that the theory is invariant with respect to conformal 
transformations of the spectral variable x. The conformal invariance in the x-plane leads to a stronger 
condition, 

dz\ 2 

— J ( T{z) ) s = analytic function of x (3.17) 

with possible singularities at the punctures x a = x{a a ). The condition (3.17) is satisfied everywhere 
on the Riemann surface except at the branch points 7 where d z x = 0. This means that the gaussian 
field is not a good approximation near the branch points. 

If the vicinity of the branch point a s one can define an involution z f-> z, such that z ^ z and 
x(z) = x(z). Up to quadratic terms in z — a s , z — a s = a s — z. In general, the involution z o z 
depends on the branch point; it is globally defined only if the spectral curve is hyperelliptic. The 
Hilbert space associated with base point z = a s splits into odd and even sectors with respect to the 
involution z -<->• z. If the Riemann surface £ is viewed as a branched cover of the spectral plane x, the 
involution zffiis the monodromy around the branch point x s = x(a s ). This is why the odd sector 
will be called twisted sector. 8 

Each sector is characterized by its stress-energy tensor. The stress-energy tensor for the even 
sector automatically sarisfies (3.17) near the branch point, while the stress-energy tensor for the odd 
sector does not. We define the twisted component of the collective field as the projection 

drf - Hz) for z in the vicinity of Gs . (318) 

The two-point function and the stress-energy tensor of the twisted field are 9 

(d$S(z) dsgcO) = £ [S W) = \B{z,2>) - \B{z,z% (3.19) 



T^(z)dz 2 



— lim 

2 z'— yz 



d^ s \z)d^ s \z') - B^ s \z,z') 



(3.20) 



The condition (3.17) represents a non-trivial constraint for the twisted stress-energy tensor 



7 Strictly speaking, we should use the term ramification points instead of branch points. According to the standard 
terminology, the points a 3 £ E where dx(z) — are ramification points while the points x(a s ) £ C are branch points. 
Each branch point is the image of one or more ramification points in the a;-plane. 

8 The richer structure of the Hilbert space at the branch points is related to the extended symmetry. For a branch point 
this is the u(2) = su(2) © u(l) current algebra. In the case of a branch point of order m the symmetry is u(m). The 
Cartan subalgebra of u(m) is spanned by the m currents associated with the different sheets of the Riemann surface near 
the branch point. 

'Since the function is symmetric with respect to exchanging the two arguments, one can antisymmetrize only with 
respect of the second argument. 
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We assume that all branch points are simple. Moreover, we assume that these points are neither 
poles of y neither zeroes of dy. If a s is one of the branch points, i.e. dx/dz ~ (z — a s ), (3.17) implies 
that the Laurent expansion of the stress-energy tensor contains only even powers of z — a s , 

TH(z) = ^41(z-a s )- 2 "- 2 , (3.21) 

and the non-vanishing Virasoro operators L2 n satisfy 

(41) =0> n>-l. (3.22) 



3.4 Dressing operators 

Our aim is to construct dressing operators located at the branch points, which restore the conformal 
invariance. For that we will apply the operator formulation (3. 15)— (3. 16) to the special case when the 
extra punctures are placed at the branch points z = a s (s = 1, . . . , n B ). 

Given a local coordinate £u s ), such that C(a B ,s) = 0, the mode expansion for the twisted field 
(3.18) is given by the odd part of the mode expansion (3.6), 



p odd 



= i E 4* 1 (3 - 23) 



h 

p>3, odd 



We denote by |0|J) the twisted Fock vacuum associated with the branch point a s . From the point 
of view of a CFT on the spectral plane, this state can be considered as the result of the insertion at 
the point x s of a twist operator of conformal weight jg. The Hilbert space space associated with the 
ensemble of the branch points of the Riemann surface is spanned by the states of the form 

(0twin4! i] > n j[ sjotw) (Pi>i, odd), (3.24) 

i i 

where the left and right Fock vacua are defined as 

(Otwl = ®<0H|, |0 tw ) = (g) |0S). (3.25) 

s=l s=l 

Once we have an operator representation of the expectation value in (3.17), the construction of the 
dressing operators becomes a purely algebraic problem. The result does not depend on the choice of 
parametrization, but the mode expansion (3.23) does. Let us pick a canonical parametrization in the 
vicinity of each branch point. We would like to choose the local parametrization variable & z a ) so tnat 

the moments //jf' are simply expressed in terms of the moduli M p s ^ of the spectral curve, defined as 10 
Mki=i yWf) , = 3,5,.... (3.26) 



"These moduli are linear combinations of the ACKM moments [11]. 
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The most convenient choice of £(z,s)> which we adopt in the following, is 



^ >s) = V2 [JWJ 3 1] V3 ^x(z)-x(a s ). (3.27) 
Then the moments (3.16) at the branch point a s are given by 

$ = {Ml s] y P/3 MW. (3.28) 
The left state (E| in the operator representation (3.14) is 

(S| = (OtwIexp^^^^Ufl + i^^HjH), (3.29) 

where s, s' take values 1, . . . , n B and p, q are positive odd integers. By convention fi p s ' = if p < 1. 
We are looking for a right state of the form 

\n) = n\o m ), n = (g>rt s \ (3.30) 

s=l 

where fiM i s th e dressing operator associated with the s-th branch point. The operator fiM is deter- 
mined from the conformal Ward identity 

(E|41^ [S] |0H)=0 (n>-l; s = l,...,n B ), (3.31) 

\s] 

where L 2 „ are Virasoro operators associated with the expansion of the twisted stress-energy tensor 
(3.20): 

T t W(z) = ^41e ( - 2 rr 4 - (3-32) 



/ft 



The explicit expression for the Virasoro operators is 

4 s i = \ E : 4 S] 4 S] ■ +^n, (n>-i). 

p+<j=2n 



(3.33) 



In order to satisfy (3.31) it is sufficient to solve the operator equation 

nH |oW) = (n > -1, s = 1,. . . ,n B ). (3.34) 



j_ 7 m 
e3h- 



Here we retained only the factor exp(^/Xg jW), with /ig = 1, from the operator deforming the left 

\s] is] 

vacuum in (3.29). This operator, when commuted to the right, shifts J p — > J p + 1^+3,0, which 
corresponds to the 'minimal' classical solution having a branch point at x s : ~ h~ l (x — x s ) 3 / 2 . 

Assuming that the dressing operator J)M can be expanded as a formal series in the creation oper- 
ators j|fp (p > 1), the solution of (3.34) is given by 



\n>l pi,...,p n 



(<?) 

S^fejW J-J J ( 3.35) 

P\...p n 11 In 
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where is a numerical factor depending on the moduli at the point a s . The sum goes over odd 
positive odd integers pj. The coefficients Wp 9 ^... tPn are universal rational numbers and are proportional 
to the correlation functions in the Kontsevich model [36]. They are nonzero only if the genus g defined 
by 

1 n 

-^ Pi = 2g-2 + n (3.36) 

6 i=i 

is a positive integer. For the sake of self-consistency, we derive in Appendix A the recurrence equation 

The numerical factor CW is fixed by the scale invariance. Upon a rescaling x — > px, the twisted 
vacuum acquires a factor p 1 / 16 and — > p~ 3 / 2 M% , while the dimensionless moments /ijf' do not 
change. The dressed twisted vacuum remains scale invariant if 

CH = (Mj s, )- 1/24 . (3.37) 

3.5 Summary: a universal formula for the h expansion 

Partition function 

The free energy Tn = In Zfj is a sum of three terms, 

J~h = F c \ + •>\ gauss + -Tint) (3.38) 

where F c \ = h^ 2 F^ is the classical action of the gaussian field on the Riemann surface, 

^gauss = J* X) = ~ In det 3 - — In M [ * ] (3.39) 

s=l 

is the sum of the gaussian fluctuations and the contributions from the scaling factors associated with 
the branch points, and F\ at is the interaction part, which vanishes when h — > 0. 
Assuming that F- iat is given by an asymptotic series of the form 

^int = J> 2 ^ 2 ^), (3.40) 

the exponent Z mt = e^' 111 can be evaluated as the scalar product 

2fat = (E|fi>. (3.41) 

The states (S | and |0) belong to the completion of the Fock space associated with the n B branch 
points and defined by the relations 

= P<W)<W, (s = l,...,n B ; pG2Z + l) 
(0 tw |jl s i = 4 Sl l°tw)=° (s = l,...,n B ; p = l,3,5,...). (3.42) 
The left state in (3.41), 



(SI d =' (0„|exp ( i £ £ + s££>W I • 0.43) 
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depends on the Bergmann kernel B and the classical field $ c i through 

B pf = j&Jto^fy^BM, p, q = 1,3,5,...; 
^ = Res y{z)dx(z) p = 5,7,.... (3.44) 

The right state is in (3.41), 

n B ( ( i\n (fl) \ 

|n> = f H«P E fi^fe^...^ , (3.45) 

8=1 \n>l ' Pl,...,p n ^l---^n y 

represents a product of dressing operators, associated with the branch points. The coefficients Wp^... tPn 
are universal numbers and are given by the genus-g, ?7,-loop correlation functions in the Kontsevich 
model (see Appendix A). 

Correlation functions 

The Fock space representation extends to the correlation functions of the collective field. They 
are given by 

< d$ qu ( Zl ) . . . d® qu (z n ) ) = < d$ qu ( Zl ) . . . d$ qu (z n ) ) gauss + < d$ qu (zi) . . . d® qu (z n ) > int (3.46) 
where ( ) gauss is the result of all possible gaussian contractions with the Bergman kernel and 

(^i)-^WL=E E sN(^)...sW(^)(S|jH...jW|n). (3.47) 

Sj=l p{>l, odd 

3.6 Diagram technique 

Performing the gaussian contractions we obtain an expression of the partition function as a series of 
Feynman graphs made of the following elements: 



Propagator: 



Tadpole: 



Vertices: 




(3.48) 



(3.49) 



(3.50) 
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The tadpoles and the vertices are associated with a given branch point and carry label s, but their 
weights are given by universal numbers. The Feynman graphs are made out of these elements by 
connecting the open lines respecting the orientation. The orientation of the lines is such that vertices 
connect to tadpoles or propagators, but not directly with other vertices. The genus g free energy is 
a sum of all connected Feynman graphs of genus g. The factors (— l) n in the dressing operator are 
taken into account by taking the tadpole with minus sign. 

The diagram expansion for the correlation functions (3.46) is obtained by adding to the Feynman 
rules the propagators 

=B(z,z>) (3.51) 



— £-<0 z =BW(z) (3.52) 

The Feynman rules (3.48)-(3.50) generalize those obtained in the one-matrix model [15] and for 
the scaling limit of various matrix models in [37, 38, 39, 40]. The structure of the genus expansion 
is such that in the case of several branch points the partition function decomposes to a product of 
Kontsevich r-functions, related by gaussian correlations. Similar decomposition formulas for the 
case of the hermitian matrix model were suggested in [41, 42, 43, 44]. 

It is possible, as is the case for the 0{n) matrix model [45], that the spectral curve has a symmetry 
relating different sheets of the Riemann surface. If such a symmetry is present, one should identify 
the twisted vacua associated with the branch points which belong to the same orbit of the symmetry. 



4 Topological recursion 
4.1 Summary 

Here we summarize the definitions and the rules of the method called topological recursion [16, 17, 
18, 19,20,21,22, 23,24, 25,26]. 

Loop observables 

Given the spectral curve (2.1), we define the n-point, genus g, functions Wn\z\, • • • , z n ) as 
follows: 

Wi°\z) = y(z)dx(z), 
wP(z 1 ,z 2 ) EE B( Zl ,Z 2 ) (4.1) 

and recursively W^^Zq, zi) = wjf^zo, z±, . . . , z n ) as 



W&fo, zj) = -£ R« K(z Q , (*, zj)W^% z IV )+ Wi 9 ~ 2 l \z, z, zi) 



■ (4.2) 



=1 m,j Jcl 



The sum inside the brackets goes over < m < g, < j < n, (m,j) / (0,0), (g,n), and all 
possible ways to split the variables zj = {zi, . . . , z n } into two sets, zj and Znj. The r.h.s. is a sum 
of terms associated with the branch points a s of the Riemann surface S. Near each branch point z ^ z 
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is defined so that x(z) = x(z). The recursion kernel K is denned in the vicinity of the branch point 

z = a s as 

*) = K(z , z) = 2{y{z) _ yCz))dx{z y * close to a s . (4.3) 
The result of the integration does not depend on the choice of the local coordinate near a s . 
Free energy 

The free energy F is defined as the formal (generically asymptotic) series 

oo 

F = Y,h 2 ~ 29 'F {9 \ (4-4) 

9=0 

whose coefficients are expressed in terms of the loop observables and the elements of the spectral 
curve, and F^ are respectively the classical action and the gaussian fluctuations of the bosonic 
field on the Riemann surface and for g > 1: 

= V Res S{z)w\ 9 \z) (4.5) 

2 — 2g Z — ' at-vn. 



z^a B 

s 



where ^(z) = y{z)dx{z). 
Homogeneity 

Under a rescaling of the one form ydx, the free energies turn into 

[Xydx] = A 2 ' 29 ^ [ydx] 

for any A G C. 
Symplectic invariance 

For (7 > 1, the free energy is invariant under canonical transformations preserving the Poisson 
bracket (2.6). 

4.2 Topological recursion in the vicinity of a branch point 

The recursion equation (4.2) involves a sum over the branch points of the Riemann surface E. If we 
take the point zq close to one of the branch points and take an appropriate scaling limit, the recursion 
kernel will consist of a single term, associated with this branch point. The piece of the algebraic 
curve which survives in the scaling limit is a curve with only one simple branch point and genus zero. 
This is the complex curve for the Kontsevich model. We will see that the solution of the recursion 
equation for the Kontsevich model plays the role of a building block for the topological recursion on 
an arbitrary spectral curve. 

4.2.1 Solution of the recursion equation for the Kontsevich model 

The spectral curve for the Kontsevich integral is a genus spectral curve with one simple branch point 
at x = xq and singular point at x = oo. It can thus be described by a rational parameterization 



x(z) = \z 2 + Xq 

y(z) = fi 3 z + fi 5 z s + fi 7 z 5 + 



< ..Liz 1 ..: a s , (4-6) 
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The moduli of the spectral curve are the moments fj,^,^,... and the position xq of the branch point. 11 
When we deform the Riemann surface, both the moments and the position of xq are varied. There is 



only one branch point at x = xq and the variable £(z,o) = ~ ^ x o = z gives a global parametriza- 
tion of the curve. 

The initial data for the topological recursion are thus the one-form 

W[°\z) = y(z)dx{z) = Hpz p ~ l dz (4.7) 

p>3, odd 

(the sum goes over p odd) and the two-form 

w(o)/ \ dz\dz 2 

W 2 y '(zi, z 2 ) = -3. (4.8) 

[zi - z 2 ) 

Let us write down the topological recursion in these terms. The recursion equation for the spectral 
curve (4.6) reads 



wQ^zk) = ~ ResK(z ,z)[ W$l(z,zj)wj££l{z,z K/ j) + (z, z, z K ) 



(4.9) 



m,j,J 

where z = —z and the recursion kernel (4.3) given by the series 

dzo 2z 



4y(z)dx{z) z 2 -z 2 
dz 1 1 



l^zlzdz l-z 2 /z 2 1 + Y (fi p /fi 3 ) Z P- s ' 

p>3, odd 

Since the residue is taken when z — > at fixed zq, we have to expand in the positive powers in z. 
Apart of the overall power l/z, the expansion is a Taylor series in z 2 . If the r.h.s. of the recursion 
equation contains only finite number of negative powers of z, which we will see to be the case, then 
the residue picks only finite number of terms. 

Let us demonstrate how the recursion equation works in the lower orders. We will evaluate 
and Expanding 

„ . , dz dz\ z z 2 

B(z,z 1 ) = ^1 + 2- + 3- + ..., (4.11) 

zf zi zf 

we find 

^3 (0) (zo,^ 3 ) = -Res K(z ,z) [B(z,z 1 )B(-z,z 2 ) + B(-z,z 1 )B(z,z 2 )] 
1 dz§dz\dz 2 



Z Q Zf Z| 



PFi (1) Oo) = - Res K(z ,z)B(z,-z) 

z— >0 

8 f4 z o ' 



(4.12) 



(4.13) 



11 Of course, the point xo can be placed at the origin by a global conformal transformation. 
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It is clear from the form of the kernel that the correlation functions with negative Euler characteristics, 
2 — 2g — n < 0, can be written under the form 



w^,...^ )= ,i- 2 '-- e n^Ep En (-£)»". ■*> -' <4 - 14) 

Pi >l,oddi=l Z i l>0 ' kj>3j=l V 3 7 

where the summation goes over odd integer p\, . . . ,p n and k±, . . . , fy. The symmetry of the spectral 
curve z — >• z/p, p2k — > P 2k+1 f^2k means that the observables (4.39) and the free energy do not change 
after this rescaling, which gives the restriction 

l n 

£(fcj-3) + l>i-3) = 3(20-2). (4.15) 

i=i t=i 

With these notations, the two examples considered above give: 

w m - 1 w {1) - - w {1) - - (4 16) 

^1,1,1 - 1 ' W 3\3 ~ g ' W 5\l - g- ^ A0) 

4.2.2 Open-close duality 

The form of the topological recursion obtained in the previous section involves both "open" moduli Z{, 
which we can visualize as boundaries on some world sheet, as well as the "closed" moduli p p (p > 3). 
The two types of moduli are related by the "loop insertion operator" [11], which depends non-linearly 
on the moments p p and the position of the branch point xq- To make our derivation self-consistent, 
we prove here that the coefficients depend on the two groups of indices in a symmetric way: 

W ki,...,k l \pi,...,p m ~ W ki,...,k k ,pi,...p m ' (4-17) 

By analogy with the string states we call this property open-close duality. 

The open-close duality follows from the study of the variations of the correlation functions when 
the initial condition ydx is perturbed. Consider the variation 

ydx(z) -»• ydx(z) + e EihlA (4. 18) 

ciz* 

for some z* G C and e small enough. This can be translated into a shift in the moduli fi^, A*5, • • • : 

fJ-p -> fJ-p - ep z~ p ~ l , p = 3,5, .... (4.19) 

Therefore the variation represents another background with the same position of the branch point and 
the same Bergman kernel, for which we can write the same loop equations. 

On the other hand, using that by definition of B = , it is easy to show [22] that the correlation 
functions transform as 

W#) (Z1) • • • , Z n ) -+ WPfa, . . . , Zn) + e W nM^--- Zn) + 0( £ 2 ). 

Comparing the two ways to write the linear term in e in the r.h.s. of the recurrence equation, we prove 
(4.17). For example, instead of (4.16), we will write 

(0) n (i) 3 (i) 5 

<1,1 = 1 » w 3,3 = o » <1 = o • 
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4.2.3 The topological recursion for the open coefficients: a trivalent theory 

Consider the simplest curve of the class (4.6), for which ^3 / ans ^5 = ^7 = • • • = 0. In this case 
the recursion kernel is 



Assume that the loop amplitudes have the expansion 

{pi >1, odd} i=l Z i 

First we notice that if all functions are even in their arguments, then only the even part of the amplitude 
W^ 0) (the Berg man kernel) will contribute to the residue. This is so if one of the arguments of 
is an external variable. If one of the arguments, say z\, is external, then we have to expand assuming 
\z\ < \z\ \ and neglect the odd part of the expansion: 

W W M= *?** _> J-jr pz P z -P, \z\<\ Zl \, (4.22) 

{z- Zl y 

where the sum goes over the odd p. Thus we can write for an expansion of the form (4.21), but 
containing both positive and negative odd "momenta" p and define the Bergmann kernel in "momen- 
tum space" as 

»S, = HW,o (p,PieZ). (4.23) 

When n = and g = 1, the first term of the recurrence equation (4.9) vanishes, while the second 
term is W^iz, —z) = B(z, —z): 

W 2 0) (Z>-Z) = ^2- ( 4 - 24 ) 

By expanding the topological recursion relation (4.9) and extracting the coefficient of ^~ 2g ~ n , 
one gets a recursion relation for the coefficients Wp) pl) ... yPn with 2g + n > 2: 



w 



q+ k= p -3 \ h=0 Jd / 

Here we used the abbreviation pj = {p\, . . . ,p n } for the set of indices I = {1, . . . , n} and the sum 
goes over all possible ways to split the set / into two non-overlapping subsets J and I\J. The solution 
of equation (4.25) can be expressed as a sum of connected Feynman-like diagrams with trivalent 
vertices. The vertex whose three legs are labelled by momenta p, q, k pointing inwards imposes the 
restriction p + q + k = 3. All internal vertices and propagators have weight 1 except the external lines, 
which have weight \p\. It will be convenient to indicate only the absolute value of the momenta and 
indicate the signs by arrows on the propagators. Then, representing the amplitudes Wp^,„ iPn as blobs 
with g handles and n legs, 
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(4.26) 



the recurrence equation takes the form 




(4.27) 



The vertex (with the given orientation of the lines) is 




(4.28) 



The black square indicates the momentum whose absolute value is put in the denominator. The prop- 
agator is the genus zero two-point function (4.23) 

rm P 1 P P 

wfl = +<>^- = \p\ 5 p+qfi = — ^K>^. (4.29) 

The initial conditions are given by the genus zero two-point function (4.29) and the genus one 
tadpole, which is computed by (4.24): 

w m= = ^o^Q) =1^3. ( 4. 30) 

Examples: 
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The recurrence equation resembles the mean field equation for a theory with cubic potential. This 
analogy can be made exact if one thinks of the string field theory Hamiltonian, which describes the 
elementary processes of splitting and joining of strings. Therefore the solution is not a series of 
Feynman diagrams, but have a tree-like structure, even though they have loops. This means that any 
diagram has a root (the starting point of the iteration) and the vertices are labeled by their height 
(the distance from the root of the tree). Because of this structure, the number of graph grows much 
slower that the number of Feynman graphs as the number of loops and external legs increases. These 
recurrence equations are of course equivalent to those obtained in topological gravity [46, 47, 48]. 



4.2.4 Topological recursion for generic values of the moduli 

Once we know the coefficients w^,„ lPn , from the general form (4.14) of the solution and the open- 
closed duality it follows that we can write down the free energy and the correlators for any background, 
characterized by the moments ^5,^7, ... by allowing lines with tadpoles. Consider the generating 
function 



w 



(g) ^Pl_ h>n 



Pl, — ,Pn 



Pi 



Pn 



n>l pi,...,p n 

The loop amplitudes (4.39) are obtained by applying the "loop insertion operator" 

d_ 



w(z 



) = Y—p- 



p>i 



to the free energy F^) [{t p }]: 

Wi 9 \ Zl ,...,z n ) 



dz\... dz n w(z 1 ) . . . w(z n )F^ 9 > [{t p - fip}] 
dz\ dz n 



p<>i 



pi+i 
1 



z 



Pn + l ' 



The new coefficients are graphically represented as follows: 




(4.31) 



(4.32) 



(4.33) 



(4.34) 



where the black blobs represent the moments /i p : 

p Hp 
« — • — 

P 



(p = 5,7,...). 



(4.35) 



The above result can be obtained directly from the recursion equation with the general recursion 
kernel (4.10). By expanding the last factor in (4.10) we obtain the recurrence equation (4.27) with a 
modified vertex, which can be graphically represented as a sum of all possible insertions: 
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(4.36) 

The same effect will be obtained if we consider each of the insertions in the dressed vertex (4.36) as 
the result of a recursion with the original vertex, with one of the amplitudes given by the tadpole 

-^0^» =- lip (4.37) 
In particular the genus g free energy for an arbitrary background is 



w 



GO 



n>l 



Pl, — ,Pn 



fan 
Pi 



Pn 



(4.38) 



4.3 Topological recursion for an arbitrary spectral curve 

Now let us consider an arbitrary spectral curve with n B branch points {a s } i=1 n as defined in sec- 
tion 2. We can cover the Riemann surface by an atlas of chartes associated with the local coordinates 
s = I, . . . ,n B . The coordinate £( 2|S ) parametrizes the vicinity of the branch point z = a s . 
In order to write down the topological recursion (4.2) in a closed form, we need to expand the 
loop amplitudes in the local coordinates. Assume that the points z\ , ■ ■ ■ , z n are close respectively to 
a si , . . . , a Sn . Then, assuming that 2g — 2 + n > 0, the n-point amplitude is expanded as a sum over 
odd integer momenta pi,...,p n : 



W^( Zl ,...,z n ) 



Pi>i 



[si...s n ] d£( ZuSl ) ^(z n , Sn ) 

Pl.«Pn C P1 + 1 £Pn + l 



(4.39) 



(2l;Sl) 



We will represent graphically the amplitudes W^ 9 > " p ™ as 



{g)[SSl ...s ni a 
ppx—Pn " p 




(4.40) 



The momenta pi are associated with the oriented external legs and the labels [sj] are associated with 
the extremities of the legs. 

Consider first a background with ijl p = 5 P ^. Then the recurrence equation for the amplitudes is 



P, / m 




P, / [S > ] 



P n P„] 




P^ ls " ] . (4.41) 



where in the right hand side one sums over s' = 1, . . . , n B and all odd integers q, k and r. The 
propagator and the vertex are denned as 

B(z,z') 



■Is'] 



B[ 



v-i 



Res Res 



(4.42) 
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p [s] Sq 



\p\ 

The propagator consists of two terms, 



Op+q+k+3,0 ^ 



Ls"] 

- . (4.44) 



P p ~ q 

The first term is the diagonal part (4.29), which conserves the momentum: 

M— * — O— « — Is"] = Csl— » — O — S— Ls'J r r l l ^/i /ic\ 

-<? = <^+g,o^, s ' b|- (4.45) 

The second term is nonzero when p and q are positive, which is illustrated by changing the orientation 
of the second line, 

[s s'] 

M h • j m = M = (4.46) 

We used the same graphical notations as in the previous chapter; the white blobs are added to com- 
pensate the factors 1/p and l/q in the definition of the black blob. 

Finally, the genus one tadpole, which is computed by (4.24) does not depend on s: 

W (l)= = ~^>^Q} =l$p,3- (4-47) 

4.3.1 Loop insertion operator 

For a generic background }, the expansion of the correlation functions near a set of branch points 

a si , . . . , a Sn reads 



«#>(*,,..., *o| = ee e ^^teStt 1 n n )«•«> 

Pi^lm^Oifcj^odd i=l ? (zi,s;) j'=l \ J 

As in the preceding case, one can add one variable in the correlation functions by using the loop 
insertion operator consisting in the shift of the differential ydx by the Bergman kernel: 

ydx{z) -> ydx(z) + e B & z *\ (4.49) 

CLZ* 

The only difference with the preceding section is that this operator acts as a shift on the times around 
all the branch points at once. Indeed, it acts as 

V [ k ~> ~ e 7fcTT~ for a11 s = 1, . . . , n B . (4.50) 

Applying the loop insertion operator to the genus g free energy, with the insertion point close to a 
branch point a s gives the open closed duality 

(g) [ Sl ...s n ]lss> s> m ] = {g) [s Sl .., n ][s' s' m ] 

pi...pn\kki...km kp 1 ...p n \k 1 ...k m • V T - /± J 



20 



4.3.2 Topological recursion for generic background 

The open closed duality allows us to get the expression of the free energy in an arbitrary background 
characterized by the moments {/i| s '} at the branch points a s ,s = 1, . . . ,n B , once we know the 
correlation functions of the trivial background fj, p = 5 P ^. Indeed, the free energies 

^>=£^ £ t^'t::: ] fl(- fJ f]- <«*> 

n>0 'pi>3,odd Sj=l i=l \ y% / 

are obtained by attaching tadpols 

is] M 

-*-m = _V2_ (p = 5 ,7,...). (4-53) 
P p 

to the external legs. 



5 Relation between the two approaches 

Observables: 

The symplectic invariants defined by (4.1), (4.1) and (4.2) are related to the correlation functions 
of the CFT on the Riemann surface (3.47) by 

d*{fr)) caaM = ^- 2+n W^( Zl ,...,z n ). (5.1) 

2g+n>2 

Furthermore, using (3.19), the recursion kernel can be written as 

l (rf$ qu (z )$lg(z)) s 
K(z 0l z ) = - rr . (5.2) 

Recursion equation and conformal Ward identity: 

The recursion equation represents an integrated form of the Virasoro constraints. To see that, we 
write the recursion equation (4.9) as the operator identity 

2d% u (z ) + V Res ( (^ qu (z )^J(z)) E gMgjMjMgA = . (5.3) 

Inserting this expression in the expectation value (5.1) one indeed obtains (4.9). 

For zq close to the branch point a s , the two-point function has Laurent expansion 

<d$ qu (*o)4u(*)> s = ^)^, S ) +re g ular - (5 - 4) 

p>l, odd 

Therefore the residue in (5.3) is a projection to the singular at zq = a s part of the integrand. The 
recursion formula (5.3) implies that 

a ^ )+ gm w igw^ (5,) 

d^(z) d^(z) 
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where <i<J>M = d<&^) +d < &qj and means 'equal up to a regular at z = a s function'. Since d<fr^ (z) ~ 
£fz s)d£(z,s) when z — > a s , (5.5) is equivalent to the conformal Ward identity (3.17). 

Graph expansions: 

The trivalent graphs that appear in the recursion procedure and the Feynman graphs of the CFT 
approach must lead to the same result. Moreover, the Feynman rules of the CFT approach can be 
obtained by partial resummation of the trivalent graphs of the recursion procedure. 

To see that, we first notice that, depending on the incoming and the outgoing momenta, only 
one of the two terms in the propagator (4.44) can be non-zero. Any trivalent graph can be split into 
connected subgraphs in such a way that the vertices belonging to the same subgraph can be connected 
only by Kontsevich (white) propagators, while the vertices belonging to two different subgraphs can 
be connected only by black propagators. Each such subgraph is associated with given branch point a s 
and contain only Kontsevich propagators (4.45). We can reorganize the sum over trivalent graphs by 
first summing up the contributions of the connected subgraphs associated with the branch points. We 
represent the sum of the contributions of all subgraphs with given values p\ , . . . , p n of the external 
momenta by an effective vertex with n legs and genus g given by (3.36). For example, 




The partial resummation recasts the trivalent theory into an effective theory whose vertices are 
universal numbers, the amplitudes Wni,...,n k of the Kontsevich model, and whose propagator and 
tadpole are given respectively by (3.48) and (3.49) with H = 1. 

6 Conclusion 

In this paper we reviewed and compared two approaches to quasiclassical expansion associated with 
a spectral curve: the CFT approach and the topological recursion (TR). The two approaches are for- 
mulated in their most general form. The only input is the spectral curve associated with the classical 
solution. Assuming that the collective field theory is conformal invariant, we derived a general for- 
mula for the genus expansion of the free energy and the current correlation functions. 

We demonstrated that the diagram techniques obtained within the two approaches, although seem- 
ingly very different, lead to the same result for the genus expansion. More concretely, we showed that 
the CFT graph expansion is obtained as a partial resummation of the graph expansion in the TR. 

The result of this paper proves that any theory solved by the topological recursion has a CFT dual 
living on the corresponding spectral curve. On the other hand, it provides a new way to prove that a 
system can be solved by the topological recursion. Namely, if a system, or matrix model, is solved 
by the CFT approach, this implies that it is solved by the topological recursion on the same spectral 
curve, which is can be more useful for explicit computations. 

The quasiclassical expansion considered here can be applyed to systems with U (N) symmetry, 
where the eigenvalues of the matrix variables have fermionic statistics and the partition function is a r- 
function. In this case the quasiclassical expansion of the free energy contains only even powers of h = 
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l/N. Our approach can be carried on also for non-trivial exchange statistics of the eigenvalues, which 
is the case for systems with O(N) or Sp(N) symmetry and, more generally, for the /3-ensembles with 
/3 7^ 1. In such systems the l/N expansion contains also odd powers in l/N. The first orders were 
found in [49], whereas a TR formalism for hyperelliptical curves was derived in [50]. In order to 
generalize the CFT method to the case of /3-ensembles, it is sufficient to add a term with a background 
charge Q ~ j3 — 1 in the definition of the stress-energy tensor (3.16), 



This will change the right state \U) in the operator representation (3.41), while the left state (S| will 
remain the same. 

Let us make some remarks concerning the application of the CFT/TR procedure to concrete matrix 
models. The case of the one-matrix model, where the spectral curve is hyperelliptic, is considered in 
details in [15]. In the case of the 0(n) model [45] and the ADE models [38] the spectral curve has a 
symmetry relating the branch points on different sheets of the Riemann surface. This makes it similar 
to a hyperelliptic curve. If such a symmetry is present, in the definition of the states (3.43) and (3.45), 
one should identify all branch points belonging to the same orbit. 

In the case of the two-matrix model with potential tr [Vi (X) + XY + V2 (Y)] , there is a symmetry 
x y and one can construct a second Fock space representation of the partition function, based on 
the conformal invariance of the spectral plane of the matrix Y. The two representations determine 
the same partition function. This non-trivial fact was proved on the TR side in [19, 51]. The first 
Fock space representation is useful for computing the conelation functions of the matrix X, while 
the second one is useful for computing the correlation functions of the matrix Y. In this context a 
generalised topological recursion procedure was developed for computing the mixed trace correlation 
functions of the matrices X and Y, or equivalently free energy in the case when the potential is 
an arbitrary polynomial of X and Y [52]. The CFT description of the 1/JV expansion for these 
correlation functions is still missing. For that one should construct a representation of the algebra 
on the complex curve. This is an interesting problem which we would like to address in the future. 
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A Solving the Virasoro constraints near a branch point 

Here we obtain the solution to the Virasoro constraints (we have put K = 1) 



T(z)dz 2 



\ lim d$(z)d$(z') -B(z,z') 





(A.l) 



k+q=2n 



which determine the coefficients w. 



'Pl,—>Pn 



. We represent the oscillator amplitudes J p as 



-p 




(P > !)■ 



(A.2) 
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Then the Virasoro operators are represented by the differential operators 

4L_2 = ^2kt k+2 d k F + tl 

k>l 
k>l 

4L P _3 = J2 kqd k d q + 2 tfkdq (p>5, odd) (A.3) 

p—k—q=3 p+k—q=3 

which act on the partition function Z[t] = Z[ti, t$, . . . ]: 

L p _ 3 Z[t] = 0, p = 1,3,5,.... (A.4) 

We are interested in solution with t 3 / 0. Then the n = 0, or p = 3, constraint is solved by 
assuming that the partion function is of the form 

Z[t p ]=t; 1/24 e F ^ /3 \ (A.5) 
We look for a solution for the free energy F, which is a formal expansion in ii, £5, £7, . . . and 1/^3 ~ 

g>0 n>0 ' fci,...,fe n >l n 

The compatibility of this expansion with (A.5) implies that the coefficients w k 9 ^ , are nonzero only 
if 

n 

= 3(2g - 2 + n). (A.7) 

i=l 

Written for the free energy F = In Z, , the Virasoro constraints become a set of quadratic con- 
straints for the derivatives 



w p = -pd p F, w P:q = pqdpd q F. (A.8) 



We have to solve for p = 1,3,... 



^ (w k w q + w k>q ) - 2 ^ t k w q + t\ S p>1 = 0. (A.9) 

p—k—q=3 p+k—q=3 



In particular, the equation for p = 3 means that the generating function 

w„ \ - z 



p>i p>i 



is scale invariant with respect to t p — > p p t p , z — > pz. 

One can solve (A.9) recursively in the genus g by expanding 



,N=E*3 2-2 ^p 9) [W*3]. (A.11) 
9>1 
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The solution of (A.9) can be obtained as a sum of graphically represented by introducing a trivalent 
vertex which gives the conservation law of the three incoming momenta p, q, k: 

v P,q,k = —$p+q+k-3,0- (A. 12) 

This vertex satisfies the selection rule (A.7) with g = and n = 3. Then (A.9) takes the form 

= V p ,- k ,- q [lJ2 W k )w q~ h) + ) + 2Vp, k ,- q t k W^ + \v pXl t\. (A.13) 

V h=0 ) 

This is the generating function for a set of recurrence equations for the coefficients w p ^... tPn . The 
equation for the correlator with n+ 1 legs p,p\, . . . , p n is obtained by differentiating w.r.t. t Pl , . . . , t Pn 
and then taking t q = <5 9) 3. The equation has the following graphical representation: 




(A. 14) 



This graphical representation is related to that for the topological recursion (A. 14) by replacing 




(A. 15) 
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